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ABSTRACT 

With recent advances in asteroseismology it is now possible to peer into the cores of red giants, potentially 
providing a way to study processes such as nuclear burning and mixing through their imprint as sharp structural 
variations - glitches - in the stellar cores. Here we show how such core glitches can affect the oscillations 
we observe in red giants. We derive an analytical expression describing the expected frequency pattern in the 
presence of a glitch. This formulation also accounts for the coupling between acoustic and gravity waves. From 
an extensive set of canonical stellar models we find glitch-induced variation in the period spacing and inertia of 
non-radial modes during several phases of red-giant evolution. Significant changes are seen in the appearance 
of mode amplitude and frequency patterns in asteroseismic diagrams such as the power spectrum and the echelle 
diagram. Interestingly, along the red-giant branch glitch-induced variation occurs only at the luminosity bump, 
potentially providing a direct seismic indicator of stars in that particular evolution stage. Similarly, we find 
the variation at only certain post-helium-ignition evolution stages, namely, in the early phases of helium-core 
burning and at the beginning of helium-shell burning signifying the asymptotic-giant-branch bump. Based on 
our results, we note that assuming stars to be glitch-free, while they are not, can result in an incorrect estimate 
of the period spacing. We further note that including diffusion and mixing beyond classical Schwarzschild, 
could affect the characteristics of the glitches, potentially providing a way to study these physical processes. 
Subject headings: stars: evolution — stars: oscillations — stars: interiors 


1. INTRODUCTION 

The cores of red-giant stars hold the key to answering a 
number of unresolved questions about fundamental physics 
that govern stellar evolution, such as mixing process, rota¬ 
tion, and the effect of magnetic fields. It has been known 
for over a decade that red giants^ show_stochastically-driven 
oscillations like the Sun (IFrandsen et all [20021) . but only 
with recent data from space missions like CoRoT and Ke¬ 
pler, have asteroseismic investigations revealed details about 
the cores of red giants. This advance has been possible 
due to the fortunate circumstance that gravity waves (here¬ 
after, g modes) in the cores of red giants couple to acous¬ 
tic waves (hereafter, p modes) in the envelope, resulting in 
mixed modes whose information about the core pro perties 
is therefore observable at the surface dDupret etalJ 120091 
[Bedding et alj l20lot) . Recent findings include the distinc¬ 
tion between stars with inert cores from those that possess 
core burning (e.g. [Bedding et aljl201lh . the measurement of 
core rotation rates much slower than predicted by current the- 
ory of angular momentum transport (e.g. iBeck et ail 120121 : 
IMosser et al J 120 1 23 ICantiello et all 1 20141) . and an ability to 
determine the evolutionary_stages of stars with unprecedented 
precision (IMosser et al.ll2014l) . Despite these findings, the full 
potential of current asteroseismic data can only be realized 
if all aspects are understood about how the internal structure 
of stars may influence the observed oscillations. To achieve 
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this, it is necessary to explore how sharp structural variations 
inside a red giant could impact its oscillation frequencies. 
Sharp structural variations can be found in stellar interiors 
at the borders of convectively mixed regions, in regions of 
ionization of elements, or between layers that have acquired 
different chemical composition as a result of nuclear burn¬ 
ing. The signatures they imprint on the oscillation frequen¬ 
cies have already been stu died observational l y and theoreti- 


cally in white dwarfs (e.g. WinJet et al.1 1991; Brassard et alJ 

1992b, main-sequence stars (Roxb 

rgh & Vorontsov 2001; 

Miglio et all 2008b Degroote et ali 

2010) including those 

like the Sun (Monteiro et al. 2000 

Monteiro & Thompson 

2005; Cunha & Metcalfe 200'S Cunha & Brandao 2011; 


Houdek & Goug h 2007; Maz umdar et alJl2014l), and in sdB 
stars (ICharpinet et al.l 20001 0stensen et al.ll2014h . which are 


essentially the cores of previous red giants. 

In red giant stars, only the signature of the helium ioniza¬ 
tion zone has been studied dMiglio et al.ll201(f) . This signa¬ 
ture arises from the stellar envelope and affects the acous¬ 
tic modes, but simulations indicate its application as a diag¬ 
nostic tool on single stars might be limited ( Broomhall et all 
120141) . However, the effect of sharp variations occurring in the 
deeper layers near the cores of red giants has neither been in¬ 
vestigated theoretically, nor been reported from observations. 
Given the mixed character of the waves in red giants, the study 
of this phenomenon requires understanding the combined ef¬ 
fect of the sharp structural variation and of the coupling be¬ 
tween p and g modes. 

Here we present the first comprehensive study of the effect 
on the oscillation frequencies of red giants from sharp struc¬ 
tural variations located in their deeper layers. We illustrate 
the impact this may have on common asteroseismic diagrams 
and investigate where this effect might be relevant during the 
red giant evolution phase. While no assumption is made about 
the degree of the modes in the analysis presented here, all ex¬ 
amples provided are for dipole modes, because these are the 
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Fig. 1.— One solar mass red-giant models considered for a detailed anal¬ 
ysis. (a) Position in the HR diagram: model la is on the red-giant branch 
just below the luminosity bump and model lb is between helium flashes.; (b) 
and (c) show, respectively, the helium profile and buoyancy frequency in the 
inner region of model la (solid curve) and model lb (dashed curve). The 
sudden decrease of the buoyancy frequency, at r/R ~ 0.055 for model la 
and r/R ~ 0.135 for model lb. marks the lower boundary of the convective 
envelope in the corresponding model. 


most promising from the observational point of view. 

2. STRUCTURE OF THE G-MODE CAVITY 

Internal gravity waves have frequencies below the buoy¬ 
ancy (or Brunt-Vaisala) frequency and propagate only where 
there is no convection. While on the red-giant branch a star 
is powered by hydrogen burning in a shell surrounding an in- 
nert radiative helium core. The g-mode propagation cavity 
extends essentially from the stellar center to the bottom of the 
convective envelope. Once stable core-helium burning starts, 
the central part of the core becomes convective, reducing the 
size of the g-mode cavity. For massive stars the transition 
between these two phases is smooth. However, according to 
current standard ID stellar models, in lower-mass stars with a 
degenerate helium core, this tr ansition involve s a s uccession 
of off-centered hel ium flashes (Bildsten et al. 2012) (see also 
Salaris et al. 2002, for a general overview of red-giant evolu¬ 


tion). 

The propagation speed of the gravity waves depends on the 
buoyancy frequency. Consequently, variations in the buoy¬ 
ancy frequency inside the g-mode cavity may perturb the pe¬ 
riods of high-radial-order modes away from their asymptotic 
value. Sharp variations in the buoyancy frequency during 
the red-giant phase usually result from local changes in the 
chemical composition. Examples of these variations are il¬ 
lustrated in Figure [T] where we show two red-giant models 
at different evolution stages (panel a), prior to and during 
the helium-flash phase, respectively, and their corresponding 
helium abundances (panel b) and buoyancy frequencies, N 
(panel c), for the core region, where N is defined by the rela¬ 
tion. 


f 1 dlnp dlnp\ 
\ 7 i dr dr ) 


( 1 ) 


Here, r is the distance from the stellar center in a spher¬ 
ical coordinate system (r,0,<p) and g , 71 , p and p are, 
respectively, the gravitational acceleration, the first adi¬ 
abatic exponent, the pressure, and the density in the 
model. The models were c ompu ted with the evolu¬ 
tion codes AS TEC (IChristensen-Dalsgaardl l2008bl) and 
MESA (iPaxton et al.ll2013l) . respectively. Two spikes are vis¬ 
ible in the buoyancy frequencies. The spikes located at rela¬ 
tive radii of « 0.003 (model la) and ss 0.005 (model lb) re¬ 
sult from the chemical-composition variation at the hydrogen¬ 
burning shell. The spike furthest out in model la, at a rela¬ 
tive radius of « 0 . 02 , results from strong chemical gradients 
left behind by the retreating convective envelope which, dur¬ 
ing the first dredge-up, extended to the region where the gas 
had previously been processed by nuclear burning . 6 As the 
convective envelope retreats, the g-mode cavity expands to 
include the sharp variation in the chemical composition; this 
eventually disappears, when reached by the hydrogen-burning 
shell which is moving out in mass as the helium core grows. 
In the case of low-mass stars, this takes place while the star 
is still on its way up the red-giant branch, when it reaches 
the well-known luminosity bump. The bump shows itself as a 
temporary decrease in luminosity when the hydrogen-burning 
shell gets close to the sharp variation in the chemical compo¬ 
sition. As a result of the decrease in the average mean molec¬ 
ular weight in the region just above the shell, the luminosity 
of the hydrogen-burning shell decreases. This is followed by 
a return to increasing luminosity when the hydrogen-burning 
shell reaches the sharp variation . 7 (Hekker and Christensen- 
Dalsgaard, in preparation). Finally, the innermost spike in 
model lb, at a relative radius of « 0.0008, results from 
the chemical composition variation caused by a helium flash. 
Spikes in the buoyancy frequency may have yet a different ori¬ 
gin from those discussed above. In particular, they can result 
from sharp variations in chemical composition left by retreat¬ 
ing convective cores that were active either during the main 
sequence or during the helium-core-burning phase. These will 
be illustrated in section [5] where we look at sharp buoyancy 
variations along the red-giant evolution more broadly. 


6 Since the model does not include diffusion, the dredge-up should leave 
behind a discontinuity in composition. However, the numerical treatment of 
the mesh in the ASTEC calculation causes numerical diffusion which leads 
to some smoothing of the composition profile and hence broadening and low¬ 
ering of the buoyancy-frequency spike, as is evident in Fig. 1. A similar but 
less pronounced effect appears to be present in the MESA models. 

7 For stars more massive than 2.2 Mq helium burning is ignited before the 
hydrogen-burning shell reaches the discontinuity and no bump occurs on the 
red-giant branch. 
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Whether or not the spikes in the buoyancy frequency are 
sufficiently sharp to produce a significant deviation of the fre¬ 
quencies of high-radial-order g modes from their asymptotic 
value depends on how the characteristic width of the spikes 
compares with the local wavelength. A comparison of the 
two scales is illustrated in Figure |2] for the two models pre¬ 
sented in Figure [I] The eigenfunction T shown in this figure 
is related to the Lagrangian pressure perturbation (see sec- 
tion[3]for a precise definition). Clearly, the width of the inner 
spike is much larger than the local wavelength in both models. 
Hence, this spike is seen as a smooth variation by the wave 
and is well accommodated by asymptotic analysis. In con¬ 
trast, at the outermost spikes the buoyancy frequency varies 
at a scale comparable to or shorter than the local wavelength. 
We therefore may expect these features - hereafter glitches 
- to change the oscillation frequencies from their asymptotic 
value. In that case, the period spacing may also deviate from 
the fixed value predicted by the asymptotic theory dTassoull 
H980h . 

3. GLITCH EFFECT ON THE PERIOD SPACING: TOY MODEL 

In this section we illustrate the effect of a buoyancy glitch 
on the oscillation frequencies and, consequently, on the period 
spacing. To accomplish that we consider first an analytical toy 
model in which the glitch is assumed to be infinitely narrow 
and well modeled by a Dirac delta function. In the analysis 
we first introduce the analytical description of the problem, 
then consider the effect of the glitch on pure g modes and, 
finally, consider the same effect when the latter couple to the 
envelope p modes. 

3.1. Setting the problem 

Our starting point for the analytic analysis is a second- 
order differential equation for the radial dependent part of 
the Lagrangian pressure perturbation, Sp, derived from the 
equations that describe linear, adiabatic perturbations to a 
spherically symmetric star, under the Cowling approximation 
(hence neglecting the Eulerian perturbation to the gravita¬ 
tional potential). This equation can be written in the stan¬ 
dard wave-equation form ( Gough|[l993l . l2007h by adopting 
\p = (r 3 /gpf) 1 ^ 2 5p as the dependent variable, where / 
is a function of frequency and of the equilibrium structure 
the f-mode discriminant defined by equation (35) of lGoughl 
20071) ). In terms of this variable, the wave equation takes the 
form, 

+ K 2 d> = 0, (2) 

with the radial wavenumber K defined by. 



Here, L 2 = 1(1 + 1) and l is the angular degree of the mode, 
c is the sound speed, and uj c and A f are generalizations of 
the usual critical acoustic frequency and buoyancy frequency, 
respectively, which account for all terms resulting from the 
spherical geometry of the problem. The exact forms of these 
quantities can be found in equations (5.4.8) and (5.4.9) of 
iGoughl (119931) . which are reproduced in Appendix B of this 
paper. The radii where K 2 = 0 define the turning points of 
the modes. These separate the regions where waves can prop¬ 
agate (where K 2 > 0) from where they are evanescent (where 
K 2 < 0). 


For typical red giants, including the models discussed in 
section [2] there are two separate propagation regions de¬ 
bited by four turning points. We denote these points as r%, 
7 * 2 , 7 ' 3 , and 7 - 4 . We note that for the models under consid¬ 
eration, ri and rq are essentially at the center of the star 
( r = 0) and at the stellar photosphere (r = R ), respec¬ 
tively. The propagation regions and the turning points r 2 , 
73 , and rq are illustrated in Figure [3] for a representative 
dipole mode in our model la, where we show the corre¬ 
sponding mode eigenfunction derived with the pulsation code 
ADIPLS (IChristensen-Dalsgaardl2008al) . The evanescent re¬ 
gion, which is located between the turning points 73 and 73 , 
separates the two propagation cavities. To its left is the so- 
called g-mode cavity while to its right we have the p-mode 
cavity. 

In practice, N (equation 0 ) is a very good approximation 
of A f everywhere except very close to the center of the star 
and in the evanescent region between the two cavities, where 
the latter diverges. These differences will be fully accounted 
for in future work (Cunha et ah, in preparation). Nevertheless, 
in the toy model presented here, we will approximate AT by 
N from the outset. Despite this and other approximations that 
will follow, our toy model retains all important features seen 
in the full numerical solutions obtained with ADIPLS and, 
as will become clear in section 0] will be important for the 
correct interpretation of the results of the latter. 

3.2. Effect on pure g modes 
3.2.1. The eigenvalue condition 

To understand the impact of a glitch on the oscillation fre¬ 
quencies it is convenient to start by analyzing a simpler prob¬ 
lem in which we ignore any coupling between the g and p 
modes. This coupling will be considered in section 13731 

In order to find the oscillation frequencies for the pure g 
modes we need to impose adequate boundary conditions to 
the solution of equation ( 0 . Towards the center of the star 
this condition is that T decreases exponentially as r goes to 
zero. Moreover, because we are ignoring any coupling and 
because the g-mode cavity is located at such depth that the 
stellar atmosphere hardly influences the solutions, the condi¬ 
tion towards the envelope also needs to be that T decreases 
exponentially for r rq. From the asymptotic analysis of 
equation 0, which ignores the effect of the glitch, we know 
that the solution, 'hj n , sat isfying the inner boundary condition 
has the form dGoughll 19931) , 

^in - T'in Kq 1/2 sin (^j K 0 dr + ^ j , (4) 

in the region ri < r <C 7 - 2 , where, following the notation of 
Gough, we have used the symbol ~ to indicate that the two 
sides of the equation are asymptotically equal. Here, ’Lin is a 
constant and the subscript 0 on K indicates that we are not ac¬ 
counting for the glitch. This inner solution is illustrated in the 
inset in Figure [3] by the continuous yellow curve. Likewise, 
the asymptotic solution to equation 0, ’L out , th at sa tisfies the 
outer boundary condition can be written as dGoughll 19931) , 

^out ~ ^outA7 1/2 sin (/ K 0 dr + ^ , (5) 

in the region ri <C r <C rq, where 'L ou t is also a constant. 
This outer solution is illustrated in Figure [3] (inset) by the 
dashed yellow curve. 
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Fig. 2.— Asymptotic eigenfunction (solid curve) and the buoyancy frequency (dashed curve) for: (a) model la and (b) model lb. The eigenfunctions have 
arbitrary amplitude and are for characteristic eigenfrequencies of these models. The arrows mark the positions of the buoyancy frequency spikes discussed in the 
text and seen also in Figure QJ. 



Fig. 3.— Normalized eigenfunction, as function of relative radius, 
for the dipole mode with frequency v = 51.20 //Hz, computed with 
the pulsation code ADIPLS for our model la. The chosen eigenfunction 
(r^ gpf) 1 / 2 ^! = r^Sp has the dimensions of energy and is normalized to 
be 1 at its maximum value. The vertical, blue dashed lines show the position 
of T 2 and 7 * 3 , the two turning points bounding the evanescent region. The 
outermost turning point, r 4 , is also shown, while the innermost turning point, 
ri, is outside the plotted range. The g-mode cavity is located between the 
unseen r\ and 7*2 and the p-mode cavity is located between 7*3 and r 4 . The 
close-up shows a comparison between the numerical (in black) and analytical 
(in yellow) eigenfunctions in a particular region, well inside the g-mode cav¬ 
ity. The continuous yellow curve represents the inner solution derived from 
equation 0 . while the dashed yellow curve represents the outer solution de¬ 
rived from equation 


Since equations 0 and (0 are both valid well inside the 
g-mode cavity, they must be the same. The requirement that 
they be the same provides the eigenvalue condition (the condi¬ 
tion that determines which oscillation (eigen)frequencies are 
allowed by the above boundary conditions). In this case, the 
eigenvalue condition translates to 



where n is a positive integer. Hence, it is this condition that 
ensures the two yellow curves match (Figure [3] (inset)). The 
phase shift that these solutions show in relation to the full 


ADIPLS solution (solid, black curve) is due to their not in¬ 
cluding the coupling to the p modes. 

Next, we include the effect from a glitch in the buoyancy 
frequency. To keep the toy model simple we will initially 
assume that the glitch appears at a single position in radius, 
r = r*, well inside the g-mode cavity, such that the asymp¬ 
totic solutions (01 and 0 are still valid on either side of it (this 
assumption will be relaxed in section U.2.31 l. Accordingly, we 
represent the glitch by a Dirac delta function, S, such that the 
buoyancy frequency becomes, 

N 2 = V 0 2 [1 +A6(r- r*)] , (7) 

where A has dimensions of length and is a measure of the 
strength of the glitch, and Nq is the glitch-free buoyancy fre¬ 
quency. By imposing continuity of the solutions 8 given by 
equations 0 and 0 at r = r* we find, 


sin (/ r ? Kodr + f) „ 

^in — / \ ^out- 

sin (/;;» + !) 


( 8 ) 


Because under the approximation considered here the glitch 
is infinitely narrow, the first derivative of the solution is not 
continuous at r = r*. The condition to be imposed on the 
derivative can be found by integrating the wave equation 0 
once in a finite region of width 2e across the glitch and then 
taking the limit when e goes to zero. Accordingly, we have. 




(9) 


where now K takes the glitch into account, differing from Kq 
only at r = r*, where N differs from Nq. Well inside the 
g-mode cavity K (equation 0 , with AT replaced by N ) may 
be approximated by. 


and, thus, we find, 


K 


LN 


( 10 ) 


dlkout dThn 
dr dr 


-AI<1 (r*) f' (r*), 


( 11 ) 


8 Strictly speaking, the continuity condition is satisfied by 5p. However, 
we have verified from the numerical solutions computed with ADIPLS that 
this condition is also very closely satisfied by \P. 
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when e —> 0. 

By differentiating equations ( 0 ]) and (0 and neglecting the 
small terms resulting from the derivatives of the amplitudes, 


K, 


- 1/2 


we 


find, after substituting in equation (fill) . 


'£out-K'o /2 (r*) cos 


'hin-?fo /2 (r*) cos 


(/: 


J 4’h out ifp / ' 2 ( r+ ) sin 


K 0 dr + - ] + 


K 0 dr + - ) = 


Kodr + ■ 


( 12 ) 


Using the continuity condition 0 and the fact that /\ (J (r* 
LNo(r*)/uir+ = LTVq/wt*, equation (fl2l> becomes 


,LN, o* . 

A -sin 

r*uJ 


(/; 


-Kodr + —) m 


Kodr + — ] sin 


>»(£ 


-Kodr + £ ) • (13) 


Equation (fl3l) provides us the eigenvalue condition in the 
presence of a glitch. We note that this condition differs from 
those deriv ed following simila r principles by Brassard et ah 
(1199211 and iMiglio et al.1 (0008) for g modes in white dwarfs 
and main-sequence stars, respectively, in particular because 
we model the glitch by a Dirac delta rather than a step func¬ 
tion. 9 

To write the eigenvalue condition in a form that can be com¬ 
pared with the one derived without the glitch, we use the rela¬ 
tion 


r r * 'n- r r 2 _ pr 

/ Kodr + — = / Kodr + — — / Kodr — (14) 

J n ^ J r\ " J r* ^ 

Introducing equation (fl4b in equation (IT3l > we hnd, after some 
algebra. 


Sill 


A'ndr + — + $ ) = 0. 


(15) 


In the above, the phase $ is defined by the following system 
of equations. 


B cos<b = 1 -A^KjL 


r*w 


sm 


A n dr + — x 


C0S (-C K ° dr + f) 


B sin $ 


= A^l 

r*u> 


sm 


Kodr ■ 


(16) 

where B is a function of frequency, also defined by the system 
of equations ( IT6b . Thus, we arrive at the final form of the 
eigenvalue condition for our toy model when including the 
glitch in the buoyancy frequency, namely, 

J K 0 dr = 7r (n - 0 — (17) 

By comparing equations 0 and (fT71 > we see that the frequen¬ 
cies of pure g modes are modified by the glitch through the 
frequency dependent phase d> only. 

9 We note that the mathematical derivation of the eigenvalue condition in 
the present work differs substantially from that presented by I Brassard et al J 
fl99l and IMiglio etakl f2008i) . in that it is based on a single equation for 
the variable VI/, rather than on the equations for variables related to the radial 
displacement and the Eulerian pressure perturbation. 


3.2.2. Effect on the period spacing 

Having considered the effect of the glitch on the g-mode 
frequencies, we now turn to the impact it has on the g-mode 
period spacing, defined as the difference between the periods 
of two modes of the same degree and consecutive radial or¬ 
ders. A possible way to proceed would be to solve the eigen¬ 
value condition numerica lly (as done, e.g., bv IBrassard et alJ 
1 19921 IMiglio et al.1 120081) to find the oscillation frequencies 
and, thus, compute the period spacings. Instead, we opt for 
deriving an analytical expression that directly describes the 
period spacings as a function of the oscillation frequency in 
terms of the glitch parameters, which we find may be a use¬ 
ful path for the future comparison with the period spacings 
derived from real data. 

Under the asymptotic approximation, the period spacing for 
high-radial order g modes, AP as , is essentially constant and 
given by ( Tassoull 1980l) . 

27T 2 

AP as ~-, (18) 

w g 

where, 

f r2 LNn 


-dr. 


(19) 


To see how the period spacing is modified from the asymp¬ 
totic value in the presence of the glitch, we first re-write the 
eigenvalue condition (fTTI i as. 


2n 


<E> 


n 4 


( 20 ) 


where P = 27r/w is the oscillation period (and we recall 
that *I> is itself a function of P). In deriving the above, we 
have used the fact that well within the g-mode cavity I\ 0 ~ 
LNo/urr to approximate J£ 2 Kodr by J£ 2 LNo/u>r dr. Be¬ 
cause Kq goes to zero towards the turning points, this ap¬ 
proximation leads to a slight overestimate of the value of the 
integral. However, it allows us to derive a simple analytical 
expression for the period spac ing. 

Next, we follow [Christensen-Dalsgaardi i ( 2012 ) 10 and define 
a function G(P), by 


(aJ cr 


S(P) = gP 


<I>. 


( 21 ) 


Using expression (l20l > and the definition of Q. we can then 
write 

7T « G (P„+l) - G (Pn) W ^AP, (22) 

where A P = P „ + 1 — P n is the period spacing in the presence 
of the glitch, or, equivalently, 

AP^. (23) 

dP 

By differentiating G with respect to P and substituting in 
equation ( |23| ) we find that this period spacing is related to the 
asymptotic period spacing by 

A P,= A P as 

A P « - -4^- = -(24) 


1 u) av 
u> s du; 


1-Pg 


10 Note that the analysis of the simplified model discussed in Section 4.2 of 
that paper contains two errors that fortuitously cancel. One is the neglect of 
a singularity in the asymptotic expression (equation (1) of that paper) in the 
evanescent region. The second is a simple sign error in the analysis leading 
to equation (22) of that paper. The combined effect of the errors is that the 
equation has the correct form, and the remaining analysis is still valid. 
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The deviation of the period spacing from its asymptotic value 
is reflected in the term Fg- Its dependence on the glitch pa¬ 
rameters can be made explicit by solving the system of equa¬ 
tions m • Defining, 


LNn 


dr 


(25) 


and making [’? A'odr ss uj* we find. 


Fg = 


and 


ALNl(<4 

r+uj g B 2 \ uj 


■ cos 


w 


+ 1 - 


ALNqUJ* 
t+uj 2 


sm 


UJ 


(26) 


B 2 = 


1 - 


ALNq 

2r*w 


cos ( 2— 

UJ 


+ 

ALNq . 2 
-sm 



L 

V w 4y J 


(27) 


In the above, the dependance of the period spacing on the 
characteristics of the glitch is expressed by the parameters A 
(glitch strength) and r* (glitch position). 

The period spacing derived from expression (124) for our 
model la is illustrated in Figure [4)i (solid curve). It varies 
around the asymptotic value (horizontal dotted line), form¬ 
ing relatively narrow dips that alternate with wider, less pro¬ 
nounced humps. The narrowing of the dips with decreasing 
frequency is due to the 1 /uj dependence of the arguments 
of the sinusoidal functions present in expressions (f26l) and 

dm 11 

Because all sinusoidal functions present in the definition 
of Fq and B (equations d26t and (l27l) . respectively) can be 
written in terms of the argument 2 uj*/uj (= 7r _1 WgP), we 
expect the distance between dips to be constant in period and 
equal to 27r 2 /uj*. Thus, it provides a measure of the depth of 
the glitch in terms of the normalized buoyancy depth. 


%. -i r ^ d , 

Wg Wg J r r 


(28) 


which is analogous to the normalized acoustic depth used 
in studies of acoustic waves 12 . Down to the middle of the 
cavity (located at Wg/w g = 0.5), the deeper the glitch loca¬ 
tion, the smaller the spacing between dips. For yet deeper 
glitches (ojg/wg > 0.5), the spacing between dips increases 
again, mirroring t he separation found for a glitch positioned 
at 1 — Wg/wg (e.g. Montgomery et ali[2003i) . 

For model la, if we take r* = 0.02164P (the radius at 
which TV 2 — TVg is maximum) we find u*/uj g = 0.0057. 
Hence the glitch is very close to the the edge of the cavity 
when measured in terms of uj^/ujg (equation (f28l> ). 

Figures and |4j; show the results of moving the glich 
deeper inside the cavity, to uj*/ui g = 0.0275 and uj*/ui g = 


11 We emphasize that unlike the case of the dips caused by mode coupling, 
the glitch-induced dips are not associated with the presence of an extra mode. 
Thus, the decrease in the period spacing at the dips is fully compensated by 
its increase at the wider, less-pronounced humps. 

12 Here we adopted the notation of IMontgomerv et all 12003 1). where 
the buoyancy depth is defined as the inverse of a period, resulting in 
the sinusoidal part of the eigenfunction having approximately the form 
sin( 7 r — jJP + 7 r/ 4 ). However, we note that the term buoyancy depth is 
sometimes used for L/u: g fe.g. lMiglio et all2008l) . instead. 
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Fig. 4.— Period spacing for pure g modes in model la. The horizon¬ 
tal dotted line shows the glitch-free period spacing, AP as . (a) Results de¬ 
rived from expression 1241 for the glitch parameters estimated for model la, 
r* = 0.02164 R (or cu*/cu g = 0.0057)) and A = 1.8 x 10 ~ 3 R (see text 
for details), (b) Results from expression (£0, but fo r a s lightly deeper glitch 
at cj*/cug = 0.0275. (c) Results from expression (24), for an even deeper 
glitch, at tu*/tug = 0.4600. The inset is a close-up of the region enclosed by 
the dashed box. Diamonds show the individual modes. See text for details, 
(d) Results obtained from integrating the wave equation numerically, ignor¬ 
ing the coupling to the p mo des (solid black curve). Overplotted is the result 
derived from expression ED for the glitch parameters estimated for model 1 a 
(dotted-dashed red curve; same as solid black line in panel a) and the result 
obtained from the same expression with the glitch parameters adjusted to 
the numerical solution, namely r* = 0.02217? (or uj*/u g = 0.0052) and 
A = 1.65 X 10 — 3 i? (red, dashed curve). The latter has been shifted in 
frequency by 1 //llz (see text for details). 


0.4600, respectively. As expected, the spacing between the 
dips at fixed frequency gets smaller as the glitch is moved 
closer to the center of the cavity (see also figures 8 and 15 
of iMiglio et al.|[2008l which show a similar effect for the g- 
mode period spacings in main-sequence classical pulsators). 
We note that in producing Figures |4j) and [4j we have also 
changed A from the value used in Figure |4ji. In Figure [dji, 
A was chosen such as to maintain the value of the effective 
glitch strength A = ALNq /r* (see right hand side of condi¬ 
tion (fl3) ) unchanged. Thus, the difference in the amplitudes 
of the patterns seen in Figures [4)i and [4]i results solely from 

















































Buoyancy glitches in the cores of red giants 


7 


the difference in the location of the glitch. For Figure 0, A 
was chosen such as to reduce the effective strength by one or¬ 
der of magnitude. In this limit of small effective strength the 
period spacing shows symmetric wiggles around the asymp¬ 
totic value, instead of the alternating dips and humps seen in 
the other two cases. Interestingly, in Figure 0 we can iden¬ 
tify a modulation of the period spacing on a scale larger than 
the separation between wiggles. This modulation is more no¬ 
ticeable when the distance between adjacent modes becomes 
comparable with the distance between glitch-induced minima. 
It is simply a sampling effect, as can be confirmed through in¬ 
spection of the inset of Figure0. We note, however, that for a 
glitch positioned at u)*/w E = 0.5, the period spacing between 
two minima is exactly twice the asymptotic period spacing, 
creating a perfect sawtooth diagram without the modulation 
seen in Figure 0. The modulation introduced by the limited 
sampling depends solely on oj*/oj g , thus providing an alter¬ 
native way to measure the position of the glitch. This is im¬ 
portant, because due to the limited frequency resolution of the 
observations, it might, in some stars, be easier to detect this 
larger scale modulation than the series of glitch-induced vari¬ 
ations in the period spacing 

Since in reality the glitch is not infinitely narrow, estimating 
the parameters r, and A from a given model requires a little 
thought. To estimate r* one may consider taking either the 
center of the glitch or the position of its maximum amplitude. 
However, to estimate A we need to consider how to trans¬ 
form the glitch in the stellar model into its infinitely narrow 
counterpart while keeping the area under the glitch essentially 
unchanged. Recalling that the Dirac 6 can be defined as the 
limit, 

5= lim _L e - (T - r * )2/e2 , (29) 

£-► 0 + ty/tt 


and taking e to be the characteristic half width of the glitch we 
find, from equation 0, 


AN 2 


N 2 



(30) 


where AN 2 = N 2 — A : ' ( j is the glitch induced deviation in the 
square of the buoyancy frequency. Taking r* to be the radius 
at which AN 2 is maximum, we estimate that e = 0.5 x 10 _3 i? 
and A = 1.8 x 10 ~ 3 R, for our model la. 


3.2.3. Numerical solution for pure g modes 

In the next step we will move to a more realistic descrip¬ 
tion of the effect from a glitch on the period spacing. Figure[2] 
shows that a Dirac 5 function is not a realistic description of 
the glitch in our stellar model. In principle, the analytical 
analysis could include a more realistic function to describe 
the glitch. However, that would have increased the complexity 
of the analysis whose main purpose was to provide a simple 
understanding of the seismic impact of the glitch. To obtain 
more realistic results we therefore solve equation (|2]» numer¬ 
ically, for the case of pure g modes, by adopting N from the 
stellar structure model. By comparing the results with those 
derived analytically, we can investigate the impact of the ap¬ 
proximations made in the analytical analysis and produce re¬ 
sults that are more directly comparable with the full numerical 
solutions from ADIPLS that will be discussed in section0 

To find the numerical solutions for pure g modes we ap¬ 
proximate K 2 in equation (0 by 


K 2 = -K II- 


N“ 


(31) 


The equation is then solved using a standard fourth order 
Runge-Kutta method with adaptive step size control and the 
eigenfrequencies are found by imposing that the solutions sat¬ 
isfy the boundary conditions ’T = 0 at r = 0 and r = R. 

The results are presented in Figure [4jl (solid, black c urve). 
Comparison with the analytical results derived in section lT.2.21 
for the glitch parameters estimated for model la (Figure 0t; 
also shown as dotted-dashed, red curve in Figure 01) pro¬ 
vides a number of interesting conclusions. First, and most 
importantly, the general form of the period spacing variation 
is similar in the two cases, reemphasizing that the effect of 
the glitch is the formation of narrow dips that alternate with 
wider, much less pronounced humps. However, it is also clear 
that both the depth of the dips and their separation in fre¬ 
quency are different in the analytical and numerical results. 
To understand these differences and their potential impact on 
glitch-parameter inferences based on the analytical model we 
adjust the glitch parameters such as to match the analytical to 
the numerical results. The new analytical solution is shown 
by the red-dashed curve in Figure 4d. The solution had to 
be shifted in frequency by 1 //Hz because the approximation 
P K 0 dr s=s ui* made earlier introduces a phase shift between 
the analytical and the numerical results. In practice, this may 
be accounted for by adding a phase to the arguments of the 
sinusoidal functions in the analytical model, thus, increasing 
the number of adjustable parameters by one. 

The rematched glitch location, r*, is almost unchanged 
(shifted by only « 50% of the glitch width), while the strength 
of the glitch is about 10% smaller. The latter reflects that the 
period spacing variations have a lower amplitude in the nu¬ 
merical results. This difference in the amplitudes and, more 
notably, the fact that they vary in opposite ways with fre¬ 
quency, is a consequence of the non-negligible width of the 
glitch. Towards lower frequency the g-mode wavelength be¬ 
comes shorter. Seen by the wave, a spike in the buoyancy 
frequency therefore appears smoother (less of a glitch). As a 
result, the amplitude of the dips in the period spacing becomes 
smaller towards lower frequency. However, in the analytical 
analysis the spike is modelled as being infinitely narrow. It 
is therefore always much narrower than the local wavelength 
and, hence, no reduction of the amplitude is seen. 

A second striking difference seen in Figure 01 concerns 
the small-scale variations that are present in the numerical re¬ 
sult, but absent in the analytical curve. Using our analytical 
model (equation (l24l) ) we found that these small-scale varia¬ 
tions would originate from a glitch at the hydrogen-burning 
shell. Given that the spike in the buoyancy frequency at this 
position is not seen as a glitch by the wave (as discussed in 
section 0 we inspected the derivatives of N and found that 
they show a high level of variation at much smaller scales than 
the local wavelength. By smoothing the derivatives and re¬ 
calculating the period spacing, the small-scale variations dis¬ 
appeared. We therefore conclude that their origin is purely 
numerical and has no physical meaning. 

3.3. Coupling with the p modes 

We now consider the same problem as in section lT2l but in¬ 
clude the coupling between the g and p modes. That requires 
replacing solution 0, valid for pure g modes, by the solution 
that accounts for mode coupling. 

When we consider that waves can propagate also in the p- 
mode cavity, the asymptotic solution to equation 0 that is 
valid well within the evanescent region is no longer an expo¬ 
nentially decaying function, but rather a linear combination 
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of an exponentially decaying and an exponentially growing 
function. The solution to equation 0 that matches the re¬ 
quired linear combination has the form (iGoughll 19931) . 


Tout ~~ Tout Kq ^ sin 


K 0 dr + -+ ip) , (32) 


in the region rq «C r «C r^, where T out is a constant and p 
is a frequency dependent phase, which is uniquely defined by 
the coefficients of the linear combination mentioned above. 
Its form will be discussed below. 

Equations 0 and (l32l) provide us the eigenvalue condition 
in the presence of mode coupling and no glitch, namely, 

£ Kodr = 7t - p. (33) 

The coupling phase p can be obtained from the eigenvalue 
condit ion derived by Shibahashi (1979) (see also lUnno et al.1 
(1989)), based on an asymptotic analysis of the equations for 
the radial component of the displacement and for the Eulerian 
pressure perturbation, under the Cowling approximation. Be¬ 
cause the oscillation frequencies must be independent of the 
variable used to expres_s_the pulsation problem, the eigenvalue 
condition derived bv IShibahashi ( 1979 ) must be equivalent to 
our eigenvalue condition (l33l i. Comparing the two we find 
(see appendix A, for details), 


p « atan 


tan 


CO — CO a 


(34) 


where q is a frequency dependent coupling factor that can 
take values in the range 0 < q < 1/4, where smaller val¬ 
ues imply a weaker coupling. Moreover, w a are the oscilla¬ 
tion frequencies that would be obtained for p modes in the 
absence of coupling (the acoustic resonant frequencies) and 

w p = (/; 3 4 c Mr) is approximately twice the asymptotic 

large separation. The corresponding period spacing, derived 
as in section U.2.21 is given by. 


A P: 


AR, 


1 co 2 dp 


A P as 
l-F c 


(35) 


co g dw 


The period spacing derived from expression (1351 ) for our 
model la is shown in Figure 0. The dips associated with the 
coupling to the p modes are equally spaced in frequency and 
located at the acoustic resonant (cyclic) frequencies (w a /27r). 
At these frequencies the denominator inside the arctan of (l34l > 
goes through zero and, as a consequence, p varies rapidly 
with frequency. The large derivative in frequency of p 
therefore makes Pc large, producing the dips in the pe¬ 
riod spacing. This is in agre eme nt with the di s cussion s by 
IChristensen-Dalsgaardl (12012l ) and iMosser et al.1 (12012bl) and 
with the period spacing derived from the analysis of real data 
for red-giant stars (e.g JBeck et al.li2()TTl) . 

Next, we add t he ef fect of the glitch. Following the same 
steps as in section COl wc find that the eigenvalue condition in 
the presence of mode coupling and a glitch is given by 


r 2 (i 

Kodr = 7r ri — 


-®-p, 


(36) 


where $ is now defined by the following system of equations, 


B cos <f> = 1 — A 


LN, 


sin 




Kodr + — + p ) x 


cos (/r? K o dr + f + V 3 ) 


B sin$ = A^F 


r*co 


sin 


K 0 dr + — + p 


(37) 

We emphasize that both B and $ now depend on p. This is 
to be expected, since the effect of the glitch on the oscilla¬ 
tions depends critically on the phase of the eigenfunction at 
the depth where the glitch is located, and that phase is influ¬ 
enced by the coupling. As a consequence, the relative devi¬ 
ation of the period spacings from the asymptotic value when 
both a glitch and mode coupling are present is different from 
what would be found by simply adding the deviations gener¬ 
ated by the coupling and by the glitch separately. This fact 
can be readily seen in the period spacing derived from the 
eigenvalue condition (l36l >. which has the form 


A P 


AP a , 


AP a , 


1 - 


co 


2 r 


dT 

dcu 


dp 

dcu 


1 - F G,C 


(38) 


As before, the deviation of the period spacing from its asymp¬ 
totic value is reflected in the second term, Pg,C> present in 
the denominator of expression d38l >. Its dependence on the 
glitch parameters can be made explicit by solving the system 
of equations (l37l) . from which we obtain 


F g,c = 

u 2 dp I | ALNp 

uig dcu \ r*wB 2 


cos 


W e- 

2 — 

co 


2p 


ALNo o 

-— sin 


ALN$ f 

r+ojgB 2 \ co 

+ (i- 


cos 


r*w 

< 
2 — 
co 


ALN*co* g 

nco 2 


"g 

CO 


7T 

4 


■p 


■ 2p 


sin 


— + T + <R 
co 4 


(39) 


where, B 2 is now given by 
B 2 = 


\_aln* 


2 r+cj 


cos 


2— +2lyJ 
CO 


ALNo n 

-— sir 

r*w 


00l 7T 

— + j+P 
co 4 


■(40) 


We see that Pg,C has two terms, each marked by a set of 
curly brackets. If we assume there is no coupling, meaning 
that p is zero for all frequencies, the first term vanishes be¬ 
cause dp/d co = 0, and the second term becomes equal to 
equation (l26l) . hence reducing Pg,C to Pq as expected. If in¬ 
stead we assume there is no glitch, which translates to A = 0, 
then the second term vanishes, and the first term becomes 
equal to Pc defined in equation d35l) . reducing Pg,c to Pc; 
again as one would expect. Finally, we consider that both cou¬ 
pling and glitch are present, but we look specifically at what 
happens at the acoustic resonance where the coupling domi¬ 
nates the expression for Pg,c- Here, the frequency derivative 
of p is very large, and hence the first term is generally much 
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Fig. 5.— Period spacing derived from the ana lytic al approach for model la. The inset shows a close-up around AP aB (d otted box), (a) case with mode 
coupling but no glitch, computed from express ion <351 . (b) case with mode coupling and glitch, computed from expression CD (black) and case with glitch and 
no mode coupling, computed from expression CD (red), adopting r* = 0.0221.R and A = 1.65 X 1CT 3 .R. 


larger than the second. However, we still have the two ex¬ 
tra cos and sin 2 ‘glitch-induced’ terms within the first set of 
curly brackets compared to just Tq, (= w 2 /w g dtp/d uf). This 
shows that the mode coupling, and therefore also the dips lo¬ 
cated at the acoustic resonance frequencies, is influenced by 
the glitch. 

The period spacing obtained from expression (l38l > is illus¬ 
trated in Figure^. Comparing with Figured, we see that the 
combined effect of the glitch and the coupling on the period 
spacing is predominantly a change in the depth of the dips at 
the acoustic resonant frequencies. Whenever a dip caused by 
the glitch coincides with a dip caused by the coupling with 
the p modes, the depth of the latter is reduced. But if a hump 
produced by the glitch coincides with the dip caused by the 
coupling, the depth of the dip increases. This behaviour is 
oposite to what would be found if the the combined effect 
were simply the sum of the deviations to the asymtotic period 
spacings caused by each effect separatly. The predicted be¬ 
haviour can be understood if we recall that the extent to which 
g modes couple to a p mode depends critically on the proxim¬ 
ity of their frequencies (assuming everything else remains un¬ 
changed, which is the case here). A glitch-induced dip in the 
period spacings means the g modes are locally more densely 
packed, as compared to the asymptotic case. Thus, if the dip 
coincides with an acoustic resonant frequency the number of g 
modes coupling to the p mode is greater, resulting in a wider 
and, consequently, shallower coupling dip. However, if an 
acoustic resonant frequency coincides with a glitch-induced 
hump, the number of g modes coupling to the p mode is re¬ 
duced, resulting in a thinner, hence, deeper coupling dip. 

4. INTERPRETATION OF FULL NUMERICAL SOLUTIONS 

In this section we consider the numerical solution of the full 
pulsation equations, including the perturbation to the gravita¬ 
tional potential, for the models introduced in section[2]and in¬ 
terpret them in the light of the results found with the toy model 
analysis presented in section [3] The full numerical solutions 


were computed with the pulsation code ADIPLS. Care was 
taken to have an adequate number of mesh points with appro¬ 
priate distribution to resolve the rapidly varying eigenfunc¬ 
tions in the g-mode cavity. 

The period spacing derived for our model 1 a from the full 
numerical solutions is shown in Figure [6] (solid curve). Com¬ 
paring with Figure [5] we see that the dips associated with the 
acoustic resonant frequencies are closer in frequency in the re¬ 
sults from the full numerical solutions than in the toy model. 
This reflects the fact that the large separation computed from 
the eigenfrequencies is smaller (by about 5% in the current 
case) than the corresponding asymptotic value, in accordance 
with the results of previous studies (e.g . IStello et ail 120091 : 
iBelkacem et al.ll201 3t iMosser et alJl2013l) . Letting aside that 
difference, we see that in the full numerical solutions the 
dips associated with the acoustic resonant frequencies show 
a depth variation resembling what is seen using our toy model 
(Figure^). Comparison with the period spacing derive d nu¬ 
merically considering only the g modes (section r3.2.3l) (red, 
dashed curve in Figure [6| confirms that the glitch in the 
buoyancy frequency is the cause of the larger-scale modula¬ 
tion seen in the full solution (see inset), and, thus, that the 
combined effect of the glitch and mode coupling is to re¬ 
duce/increase the depth of the dips positioned at the acous¬ 
tic resonant frequencies when they coincide with a glitch- 
induced dip/hump. 

To illustrate how a glitch in the buoyancy frequency 
could be revealed in observational data. Figure [7] shows 
(a) 1 /x/Qn i representative of relative mode amplitude 
(tChristensen-Dalsgaardi2004llAerts et al,ll2010l) . and both (b) 
a frequency- and (c) a period-echelle diagram corresponding 
to model 1 a. Here, Q n i is a measure of the inertia of a mode 
of radial order n and degree l, relative to that of radial modes 
defined by 

Qnl = (41) 

H=0 
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Frequency (p,Hz) 


Fig. 6.— Period spacings derived from ADIPLS for model la, including 
the effects of the glitch and of the coupling between the g-modes and the p- 
modes (solid curve), compared to the integration of the wave equation o. 
ignoring the coupling with p modes (red, dashed curve). The inset shows a 
close up around AP as (dotted box). 


where I n i is the surface normalized mode inertia, 
T _ Iq ‘ [£r 2 + l(l + l) Cg] pr 2 dr 

Ul M£ r ( R s ) 2 


(42) 


and Ii —o is obtained by interpolating /„o to the frequency of 
the mode under consideration. Moreover, £ r and £h are the ra¬ 
dial and horizontal components of the displacement, respec¬ 
tively, M is the stellar mass, and R s is the surface radius. 
The frequency echelle shows the frequency spectrum (Fig¬ 
ure [7}i) divided into segments of fixed length that are stacked 
one above the other. The length of segments equals the av¬ 
erage frequency separation between overtone radial modes, 
Au, found as the slope of_a linear fit to the radial modes 
versus their order dGrec et al .11 1 9831) . In the period-echelle 
diagram we show only the dipole modes, and here the ab¬ 
scissa is the mode period modulo the asym ptotic period spac¬ 
ing, AP as (equation (fl~8l >, see e.g. JBedding et al] (1201 1 Il f. For 
clarity, we show only modes of relative amplitude above 5% 
of the radial modes in the echelle diagrams. Noise set aside, 
the result seen in this model is a broadening of the clusters of 
‘observable’ dipole modes where the location of the dips co¬ 
incides with that of a cluster. We see this effect in Figure [TJi 
near 28/yHz and 41 //Hz (see also Figure [6}. In the period- 
echelle diagram, the same effect shows as a strong distortion 
of the usual ‘S’ shaped mode pattern seen between each radial 
mode order in a glitch-free case (see for example figure 1 of 
iBedding et all ()201 it) ). 

Next we consider the full numerical solutions for our 
model lb, located in the core-helium-flash evolution phase. 
The period spacings derived from the ADIPLS results for this 
model are shown in Figure [8] (solid curve). The existence of 
closely-spaced pronounced dips in the period spacing makes 
it harder to identify the dips associated with the acoustic reso¬ 
nant frequencies in this case. To help with that identification, 
we mark the frequencies of the radial acoustic modes (verti¬ 
cal lines) and recall that the dips produced by the coupling 




1 2 3 4 0 10 20 30 40 50 60 

Frequency MOD 4.48 (/zHz) Period MOD 67.55 (s) 


Fig. 7.— Model la (iz max = 45 Hz) - (a) Pseudo amplitude spectrum 
based purely on mode inertia relative to the radial modes (1 /VQnO- (radial 
in black and dipole in red), (b) Echelle diagram including radial (circles) and 
dipole (triangles) modes. The abscissa is the frequency modulo Au. Symbol 
size follows the peak heights from panel (a). The arrows indicate clusters 
of dipole modes affected by a glitch-induced dip in the period spacing, (c) 
Period echelle following the notation of panel (b). The abscissa is the period 
modulo AP as - The solid curve connects all the dipole modes. 


between p and g dipole modes should be positioned roughly 
mid way between consecutive radial modes. Indeed, single 
or double dips of greater depth than their neighbors are found 
at the expected frequencies. Comparison of the period spac¬ 
ing derived from the full solutions (solid curve) with that de¬ 
rived considering only the g modes (red, dashed curve) shows 
that the two are similar everywhere, except at the frequen¬ 
cies of these more pronounced dips. This confirms that the 
more pronounce dips are produced by the coupling between 
the p and g modes and excludes that this coupling is the cause 
for the other dips. Using our analytical model for the case 
of having a glitch but no coupling (equation (1241 0 we find 
that the less pronounced dips are caused by the outer spike in 
the buoyancy frequency (Figure [2j)), that is, the glitch at the 
hydrogen-burning shell. Moreover, we also confirm, based on 
the analytical model, that despite the glitch being located rela¬ 
tively far from the center of the cavity (at uj*/oj g = 0.075)the 
larger-scale modulation seen in the period spacing (on a scale 
of about 10 glitch-induce d dips) is explained by the sampling 
effect discussed in section U.2.21 
The separation between glitch-induced dips in model lb is 
similar to the width of the dips associated with the acous¬ 
tic resonant frequencies, which makes it difficult to interpret 
the combined effect of the glitch and the mode coupling in 
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Frequency (p,Hz) 

Fig. 8.— Period spacings derived from ADIPLS for model lb, including 
the effects of the glitch and of the coupling between the g-modes and the p- 
modes (solid curve), compared to the integration of the wave equation o. 
ignoring the coupling with p modes (red, dashed curve). The vertical, lines 
indicate the frequencies of the radial modes. 


this case. Nevertheless, the comparison between the dips at 
w 18.1/iHz and 23.2/rHz, indicates that the combined ef¬ 
fect is the same as for model la. The latter dip is placed at 
a glitch-induced hump and, consequently, has its depth in¬ 
creased, while the former dip is placed at a glitch-induced dip 
and has its central depth reduced, forming a double-dip struc¬ 
ture. The observational impact of these double-dip structures 
will be discussed in section^ 

5. OBSERVING GLITCHES IN RED GIANTS 

Next, we search an extensive set of stellar models of various 
masses to locate the stages of evolution where one could po¬ 
tentially observe the seismic signature from buoyancy glitches 
in red giants. Our stellar models are derived using the ‘de- 
fault’ work inlist of MESA-v5271 dPaxton et alJl201 ll 1201 3h 
with the only change that we turned off mass loss. These 
canonical models do not include diffusion or extra mixing be¬ 
yond convection defined by the classical Schwarzschild crite¬ 
rion (ISchwarzschildl 19061) . Our search comprises tracks rang¬ 
ing 1.0-3.0Mq, all roughly with solar abundance, spanning 
the entire evolution from the bottom of the red-giant branch to 
near the end of the asymptotic-giant branch. To check that we 
obtain consistent results we also derived ASTEC tracks for 
I.OMq to near the tip of the red-giant branch and for 2.4M 0 
to the end of helium-core burning. The frequency calcula¬ 
tions based on the full numer ical solution shown in this sec¬ 
tion were made using GYRE (iTownsend & Teitlerll2013l> . but 
spot checks were made to verify that these results were con¬ 
sistent with what we obtained using ADIPLS. 

For each model, we calculate the dipole g-mode frequencies 
by solving equation ([2} numeric ally, wi th I\ defined by equa¬ 
tion (fTil) as described in section [3.2.31 hence, neglecting the 
coupling to the acoustic cavity. The calculation is restricted to 
within a ;z lnax /2-\v'ide range centered around the solar-scaled 
Gnaxtx g/VTeff, where T e g is the effective temperature and 
i/max is the frequency of maximum oscillations power. This 
range is roughly equal to the full width at half maximum of the 


excess power observed for s olar-like oscillations (IStello et al.1 
l2004t iKieldsen et all 120051 : iMosser et all 1201 2ah . From the 
resulting frequencies we then derive the series of pairwise pe¬ 
riod spacings, A P, and calculate an index of glitch-induced 
variation in A P to determine if the g-modes are effected by a 
glitch. We tested two indices, both showing consistent results. 
One was simply the RMS of the period spacings and the other 
was the height of the strongest peak in the Fourier transform 
of the series of period spacings versus period. The latter is 
shown for a section of the 1M 0 track in Figure [9j a) indicat¬ 
ing a region bracketed by the vertical dotted lines where the 
index is above twice the floor level. This phase is therefore 
identified as showing excess variation in A P. Following the 
approach described in relation to Figure [TJt (section 4), we 
also derive the large frequency separation for radial modes, 
Av. 

5.1. Before helium ignition 

Along the red giant branch we find glitch-induced varia¬ 
tions only at one particular phase in evolution lasting roughly 
5-10 million years. Interestingly, this coincides with the lu¬ 
minosity bump. Figured b-h) summarizes the results near the 
bump for low-mass models (1.0 < M < 2.2M 0 ). Models 
beyond 2.2M 0 do not show the bump because the glitch from 
the first dredge-up is not reached by the hydrogen-burning 
shell until after the model is past the tip of the red-giant 
branch. The thick red curves indicate the phases of excess 
variation in A P. This excess variation can be attributed to the 
glitch left by the dredge-up as illustrated by model la (Fig¬ 
ure^). The only exception to this picture is along the 2.OM 0 
track, which shows an extra slightly earlier phase of excess 
variation arising from a subtle but interesting combination of 
effects, also resulting in the extra luminosity bump we see for 
this mass at 1.026 Gyr. During the main-sequence phase the 
gradually retreating convective core leaves a steep gradient in 
molecular weight (hence a spike in the buoyancy frequency) 
where the convection reached its maximal extent at young 
age. For models below 1.8M 0 , the gradient is smoothed away 
by the hydrogen-burning shell, which is later established at 
almost the same location. However, for models of roughly 
2.OM 0 , the hydrogen-burning shell starts at a smaller radius 
relative to this gradient, and the gradient therefore survives for 
a while. This allows the star to evolve to the point where the 
local wavelength becomes comparable to the scale of the asso¬ 
ciated spike in the buoyancy frequency, giving rise to the first 
phase of excess variation in A P that ends when the hydrogen¬ 
burning shell finally reaches the location of the spike. In more 
massive stars that same spike is erased by the first dredge-up 
before the spike appears as a glitch for the gravity waves. 

5.2. After helium ignition 

In Figure [TO] we show the result for post-helium ignition 
tracks with masses l.OM 0 , 1.6M 0 , 2.2M 0 , and 2.8M 0 . 
Again, thick red curves indicate evolution phases showing 
excess variation in A P. The downward-pointing arrows in¬ 
dicate when the last off-center helium sub-flash and associ¬ 
ated convection zone reaches the center, signifying the start of 
quiescent helium-core burning in the models with degenerate 
cores before helium ignition (Figure [TOh.c). There is no such 
equivalent for the higher-mass models, in which a more gen¬ 
tle at-center helium ignition starts immediately at the tip of the 
red giant branch. The upward-pointing arrows mark the end of 
helium-core burning at the so-called asymptotic-giant-branch 
bump, and the subsequent asymptotic-giant-branch phase. In 
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Fig. 9.— Panel (a): Maximum signal in the Fourier transform of the se¬ 
ries of period spacings along the IMq track. Vertical dotted lines bracket 
the region of glitch-induced variation in A P. Panels (b-h): Close-up of the 
evolution near the red giant branch luminosity bump as a function of age 
for models with fully or partially-degenerate cores. Thick red curves indi¬ 
cate phases of exces s var iation in A P. A MESA- equivalent of model la 
discussed in sections [2j4]is shown, in panels (a) and (b). 


Fig. 10.— Early (left) and late (right) stages of helium core burning. Thick 
red curves indicate phases of excess variation in AP. The time in Myrs since 
helium ignition at the tip of the red giant branch is indicated along the top 
axis of each panel. Down/Up-ward pointing arrows show the start/end of 
quiescent helium-core burning (see text). The left- and right-side annotation 
of the ordinate applies to both panels. 


the following we will discuss each phase in turn where we see 
excess variation in AP. 

5.2.1. Low-mass stars 

Along the 1 .OM 0 track we see repeated intervals of excess 
variation during the initial helium sub-flashing phase. Each 
of these intervals are interspersed by short off-center helium 
burning sub-flashes where the g-mode cavity is split in two 


(iBildsten et al.ll2012h . If both g-mode cavities are taken into 
account, the resulting effect on AP during this cavity s plit 
differs significantly from what is presented bv IBildsten et al.l 
( 2012 ). who ignored the inner cavity in their analysis. The 
results including both cavities will be discussed in a forth¬ 
coming paper. The intervals with only one g-mode cavity 
are illustrated by model lb, discussed in sections [2]|4] and 
model 2 (Figure ITOh). In Figure QT] we show a multi-faceted 
view of model 2 including its core structure and the glitch 
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effect on the observed frequencies. Model 2 is similar to 
model lb except that it has a lower luminosity, hence larger 
z/ max and is therefore more likely to represent a case where 
AP can be measured in observational data (Mo sser et al.1 
l2014llGrosiean et al.l2014l) . As in model lb, we see a series of 
glitch-induced dips in AP (Figure fTTbl. The associated dips 
in mode inertia, or peaks in amplitude (Figure fTTt). suggest 
that some modes would be observable even if they are far from 
the acoustic resonant frequency. This decrease in the inertia 
arises because some, almost pure, g-modes are trapped in the 
outer part of the g-mode cavity. As a result, we see a split of 
the l = 1 ridge in the echelle diagram (Figure fTTlll. That split 
is most evident where one of the glitch-induced dips coincides 
with an acoustic resonant frequency (a coupling-induced dip), 
splitting the coupling dip into two (Figure fTTbl. 

In the following quiescent helium-core burning phase we 
see no significant variation in AP for our canonical models. 
However, towards the end of core burning (Figure [TOb). the 
retreating convective core leaves a sharp glitch, which results 
in very high-frequency variation in AP at the asymptotic- 
giant-branch bump and the early helium-shell burning phase. 
Figure [T2l shows the buoyancy frequency of model 3, which 
is representative for models in this phase. The glitch is lo¬ 
cated closely to the center of the cavity (at w*/w g ~ 0.5). 
Hence, the induced period-spacing variations occur over a 
scale comparable to the separation between two consecutive 
modes, which results in a low-frequency modulation in AP 
on top of the high-frequency variation, as discussed in sec¬ 
tion [3T2T2] and illustrated by the analytical result in Figure 0J;. 
The numerical result including adiabatic frequencies and in¬ 
ertias of such models is still under investigation and will be 
presented in a forthcoming paper. The AP variations indi¬ 
cated at the early helium-shell burning phase (Figure [TUb. d, 
f, h) are similar for all masses that we investigated, and orig¬ 
inate from the same physical reasons as discussed above for 
model 3. Moreover, all models that ignite helium in a degen¬ 
erate core, which include the models shown in Figure ITOb . 
show quite similar behavior to the 1 .OMq case, and will not 
be discussed further. 

5.2.2. High-mass stars 

Moving on to a case where helium ignites in a partially 
degenerate core, we see excess variation in the early stages 
of helium core burning as illustrated along the 2.2 Mq track 
(Figure ITok). This variation originates from the glitch at the 
hydrogen-burning shell. We show a representative model in 
Figure fl3l 

Finally, representative of stars igniting helium in a non¬ 
degenerate core, the 2.8 Mq track shows two phases of ex¬ 
cess variation during early stages of helium-core burning (Fig¬ 
ure [Tot?). The first phase is the ‘high-mass’ non-degenerate- 
core equivalent to what we saw in the low-mass degenerate- 
core models near the red-giant-branch bump, where the glitch 
from the first dredge-up ‘enters’ the g-mode cavity (Figures[9] 
andGi). As in the red-giant-branch bump cases, the AP vari¬ 
ation vanishes when the hydrogen-burning shell reaches and 
smooths out the glitch, but here this occurs after the model 
has become a quiescent helium-core burning clump star with 
Az/ ~ 7.5/iHz. Figures HU [15] and [T6l show three exam¬ 
ples along this phase of evolution. Like in model 2, the 
glitch in these three models is expected to cause relative high 
amplitudes in the frequency spectrum for the almost pure g 
modes located at dips in AP (panels b and c). The echelle 
diagram can therefore appear to show a dominant spacing be¬ 


tween strong modes that is significantly larger than the under¬ 
lying period spacing between adjacent modes (Figure ITStl). 
A second phase of variation occurs due to a glitch that was 
built up near the edge of the convective core during helium 
ignition and its subsequent maximal extent. However, we do 
not show an example of this phase because, in our models, the 
variation only shows up with a relatively low amplitude (AP 
~ 10 sec) and dips in AP that are rather broad and widely 
separated, making it very difficult to detect when the coupling 
to the acoustic modes is included. 

5.3. Discussion 

Due to the glitch-induced variation in AP around AP as , 
one can choose to use either AP as (horizontal dotted line 
in panel b) or the maximum period spacing to generate the 
period echelle (see for example Figure 1 1 5b). We chose to 
use AP as , which in some cases creates one overall ‘S’ shape 
per radial mode order as in the glitch-free case (see high- 
frequency end of Figure |7J see also Fig.lc of iBedding et al.l 
( 201 li) ). but modulated with the glitch-induced variation on 
top as in Figure HB and ri3l z. Had we chosen to use the 
maximum value of AP, we would obtain one ‘S’ shape for 
every glitch-induced dip in AP. In other cases using AP as 
makes the period echelle look very complicated with no clear 
pattern, such as in Figures [TTb. while the maximum value ap¬ 
pears to create a better aligned echelle, by straightening the 
zig-zag pattern. The latter might therefore be misinterpreted 
as the asymptotic period spacing of a glitch-free star. This 
could potentially explain some of the more massive stars with 
observed AP as reported to fall significantly outside the main 
ensemble in the AP as -Az/ diagram bv lMosser et ali ( 2014 ?). if 
indeed real stars share the frequency behavior shown by our 
models. 

In search for excess variation in AP, as summarized in Fig¬ 
ures [9] and [lOl we deliberately ignored the coupling with the 
envelope (p-)modes to simplify and speedup the process. In 
real data, one would have to separate the variation in AP 
caused by the coupling from the effect induced by the buoy¬ 
ancy glitches. We verified that our results presented here are 
consistent with what we obtain if the p-mode coupling is in¬ 
cluded, which we did by first fitting and removing the cou¬ 
pling pattern from the mixed-mode frequencies derived from 
the full numerical solution (e.g. ADIPLS/GYRE), and sub¬ 
sequently deriving the RMS of the residual period-spacing 
variation. Fitting and removing the coupling pattern was per¬ 
formed along I.OMq and 2.4M . tracks using the toy model 
for the coupling presented by IStellol (12012l ). but could as 
well be done using equation (l35l > (see also equation (9) in 
iMosser et al.ll2012bh . The analysis of glitch-induced varia¬ 
tions in AP from real data will be presented in Stello et al. 
(in preparation). 

Although we verified that the results summarized in Fig¬ 
ures [9] and [TO] are similar when based on ASTEC, our 
ASTEC models generally showed less high-frequency vari¬ 
ation in AP due to the glitches being smoother, arising from 
numerical diffusion in ASTEC, as described in section[2] It 
is also expected that including additional mixing processes 
could affect the buoyancy frequency significantly, and hence 
alter the signature in the frequencies, which would potentially 
provide a way to test various prescriptions of mixing (Con¬ 
stantino et al., in preparation). 

6. CONCLUSIONS 
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Fig. 11 . — Model 2 (M = 1.0 Mg ; i' max = 31 /'Hz) - (a) Buoyancy frequency with key features indicated, (b) Period spacing of pure g modes (as in 
section [ 3 . 2 . 2 ) (red) and full numerical solution using GYRE (black). The vertical dotted lines i ndicat e the appr oximate posi t ion of t he dipole acoustic mode s. 
They have been positioned relative to the nearest radial mode in agreement with Stello et al. ( 2014 ) (see also Huber et al. 1 2010 ); Montalban et al. 1 2010 )). 
Horizontal dotted line marks A P as . (c) Pseudo amplitude spectrum based purely on mode inertia, normalized to the radial modes. Dipole modes are shown in 
red and radial modes in black, (d) Echelle diagram. The abscissa is the frequency modulo Aw. Symbol size follows the peak heights in panel (c). (e) Period 
echelle diagram. Symbol sizes as in panel (d). The abscissa is the period modulo AP aa . Dotted lines indicate the approximate position of the dipole acoustic 
modes. Black curve connects all dipole modes. 



r/R 


Fig. 12.— Model 3 (M = 1.0 Mg ; /z rnax = 15 //Hz) - Buoyancy fre¬ 
quency with key features indicated. 


We have shown that structural glitches in the cores of red 
giants can significantly affect the adiabatic properties of their 
mixed modes - both mode inertias and frequencies. The 
modulation in mode inertia can have strong consequences for 
which modes are observable. Moreover, the change in the 
frequency pattern shows up as a variation in the underlying 
period spacing of pure g-modes around the fixed asymptotic 
value of the glitch-free case. Hence, assuming the period 
spacing follows the simple glitch-free asym ptotic behavior 
(equation (l35l) . see also equation (9) in iMosser et alJ(l2012bl) '). 
can hamper the estimate of the asymptotic (glitch-free) period 
spacing, AP as . This might explain some of the stars observed 
to show a period spacing that does not follow the main ensem¬ 
ble of stars both along the red-giant branch and the red clump 
(tMosser et al.ll2014l) . 


We provide an approximate analytical solution to the wave 
equation in the presence of both a structural glitch and the 
coupling between p and g modes. We find that the combined 
effect of a glitch and mode coupling is not merely the sum 
of the two. The combined effect is a modulation of the depth 
of the dips at the acoustic resonant frequencies and, in some 
cases, the split of these in two. The glitch-induced variations 
in the period spacing are equally spaced in period, and reflect 
the depth at which the glitch is located, while the amplitude 
of the variation is a measure of the effective strength of the 
glitch. 

From an extensive set of evolution tracks of varying mass 
we find glitch-induced variation at the red-giant-branch lumi¬ 
nosity bump, at the early phases of helium-core burning, and 
at the asymptotic-giant-branch bump, which signifies the be¬ 
ginning of helium shell burning. We note that some of these 
evolution stages last for a relatively short period of time, mak¬ 
ing the detection of glitches in such stars a strong indicator of 
relative age. 





































































Buoyancy glitches in the cores of red giants 


15 



Fig. 13.— Model 4 (M = 2.2 Mg ; i/ m ax= 84 ji Hz) - notation as in Figure [TT1 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 0 50 100 150 

Frequency MOD 3.23 (/zHz) Period MOD 169.6 (s) 


Fig. 14.— Model 5 (M = 2.8 Mg ; i/ nm x= 34//Hz) - notation as in Figure [TT1 


APPENDIX A 

In this appendix we derive the explicit form of the cou¬ 
pling phase, ip, that appears in equation ( |33| >. This phase is 
uniquely determined by the coefficients entering the solution 
of the wave equation in the evanescent region J '2 < r 13 . 
In principle these can be determined by matching the solu¬ 
tion in the evanescent region to that in the p-mode cavity and, 
subsequently, applying an appropriate boundary condition at 
the photosphere. However, in red-giant models such as those 
under study, we find that K 2 defined by equation (0 goes to 
—oc at some critical radius r c located in the evanescent re¬ 


gion between the two cavities. The analysis of the wave equa¬ 
tion across this singularity is rather cumbersome and will be 
considered in a separate paper (Cunha et ah, in preparation). 
Here, we use, instead, the eigenvalue condition presented by 
IShibahashi ( 1979 ). which accounts for mode coupling but not 
for rapid variations in the structure (hence no glitch). Their 
eigenvalue condition is derived through the asymptotic anal¬ 
ysis of the pulsation equations, under the Cowling approxi¬ 
mation, for two pulsation variables, one related to the radial 
component of the displacement and the other related to the 
Eulerian pressure perturbation. The simultaneous use of the 
two equations allows the author to avoid having to match the 
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solutions across critical points similar to that referred above. 
The result is the eigenvalue condition (IShibahashil 19791 equa¬ 
tion (31)) 


cot 


/to dr I tan 


n 0 dr = q 


(43) 


where q is often called the coupling factor and is given by 


9 = - exp [ -2 


c r 3 


K 0 | dr 


(44) 


where, as before, we used the subscript 0 to indicate that this 
condition is valid in the absence of a glitch. In the above, /cq is 


an approximation to the radial wavenumber appearing in the 
equations used by the author and is given by 


K 


2 _ 
o — 




(45) 


Inside the g-mode cavity kq ~ Ljr -^/(l — Nq/lu 2 ) ~ Kg. 
Using this fact, we can combine the conditions (f33l> and ( 143b 
to write. 


tan (ip) = 



(46) 
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Next, we note thatthe eigenvalue condition for pure p modes In the above, Af is the generalized buoyancy frequency 

derived bv iShibahashil 11 979 1 (his equation (26)) is given by 


K 0a = mn, 


(47) 


where m is an integer and the subscript “a” was added to in¬ 
dicate that this condition provides what would be the eigen- 
frequencies of acoustic waves in the absence of coupling. 

Writing kq = Ko a + 5 k and taking 5 k ~ 5ui/c (which is 
a good approximation throughout the p-mode cavity, except 
near the turning points r 3 and r 4 ) we then have 


tan 



s=s tan 


u - UJ a 

Wp 


(48) 


where w a are the eigenvalues that would be obtain for p modes 
in the absence of coupling and w” 1 = f t ' j 4 c _1 dr. Finally, 
using (f48b in equation (146b we hnd 


<P 


atan 


tan 


U! — 0J a 


(49) 


APPENDIX B 


Below we reproduce the expressions for the generalized 
buoyancy frequency and critical frequency that appear in 
equati ons (5 . 4.8) an d (5.4.9) derived bv lGoughld 19931) . 

In iGoughl (1 993 1. the author expresses the equations de¬ 
scribing linear, adiabatic pulsations in terms of the Lagrangian 
pressure perturbation 5p. After performing the Cowling ap¬ 
proximation, the resulting second order differential equation 
for 5p is reduced to the standard form by defining a new de¬ 
pendent variable 41 = (r 3 /gpf) x / 2 5p, where / is the f-mode 
discriminant given by 


/ = 


o 

to r 


+ 2 + 


L 2 9 


(50) 


and H g is the scale height for the gravitational acceleration 
obtained following the general definition adopted by the au¬ 
thor that the scale height for a quantity q is H q = — . The 

wave equation resulting from this variable transformation is 


d 2 'h 

dr 2 " 


KH = 0 , 


with the radial wavenumber K defined by. 



(51) 


(52) 



(53) 


where h is the scale height for g/r 2 and is related to H g by 
h~ l = H- 1 + 2r -1 and 77 is the scale height for gpf jr 3 
and is related to other relevant scale heights in the analysis 
by 77 _1 = Hp 1 + HJ 1 + h -1 + r -1 . Moreover, w c is a 
generalization of the critical frequency given by 


uj 2 = 


m 2 


1-2 


dr 


£ 

h 


(54) 


Using the definition for the density scale height and the 
equation for hydrostatic equilibrium, the buoyancy frequency 
defined by expression 0 can be written as, 


N 2 =g 




(55) 


Comparing this expression with expression (l53l > we see that 
the generalized buoyancy frequency has 77 in the place of 
H p and includes an additional term, —2 /h. As mentioned 
bv lGoughld 1993b (and seen from the definitions of h and 77), 
these differences result from the geometry and self-gravity of 
the equilibrium state and, consequently, AT 2 reduces to N 2 
in the limit of a plane-parallel envelope under constant grav¬ 
itational acceleration. Comparison of the generalized criti¬ 
cal frequency with the one derived in that same limit ( Gough 
12007b . shows that the difference between the two is also solely 
the outcome of geometry and self-gravity. 
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